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Abstract. We study the influence of the scheme for the correction for spurious center-of-mass motion on 
the fit of effective interactions for self-consistent nuclear mean-field calculations. We find that interactions 
with very simple center-of-mass correction have significantly larger surface coefficients than interactions 
for which the center-of-mass correction was calculated for the actual many-body state during the fit. The 
reason for that is that the effective interaction has to counteract the wrong trends with nucleon number 
of all simplified schemes for center-of-mass correction which puts a wrong trend with mass number into 
the effective interaction itself. The effect becomes clearly visible when looking at the deformation energy 
of largely deformed systems, e.g. superdeformed states or fission barriers of heavy nuclei. 
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1 Introduction 

It is generally known that the ground-state wave func- 
tions of mean-field models break symmetries which had 
been originally given in the many-body Hamiltonian or 
effective energy functional. Violation of translational sym- 
metry is unavoidable because the center-of-mass of a sys- 
tem is localised by the mean-field potential. This causes 
a spurious contribution from the center-of-mass vibra- 
tions to the energy and other observables. The problem 
has been discussed since decades and several solutions 
have been developed in the course of time, for an overview 
see Q]. A rigorous way to restore the broken symmetries 
is the projection method. Projection-before-variation is 
the perfect solution which has been applied even in re- 
alistic applications |^, but it constitutes a numerically 
extremely challenging task and is still too costly to be 
used in large-scale investigations of nuclear structure. A 
simpler approach is the projection-after-variation method 
where the mere HFB state is varied but projected wave 
functions are used to calculate observables W . A study of 
this approach in the context of self-consistent models has 
hinted that full projection effects could be quantitatively 
important in light nuclei |^. Nonetheless, by far the most 
applications deal with approximate ways to compute the 
center-of-mass corrections for reasons of feasibility and 
transferability. The standard procedure is to expand the 
correction in orders of moments (Pc"n ) ^^'^ to stop at first 
order W . And even that is often further simplified in var- 



ious manners. As a consequence there are several recipes 
around for performing the center-of-mass correction. 

This diversity as such could possibly be bearable. The 
situation is complicated by the fact that all quantita- 
tively successful nuclear mean-field theories employ a phe- 
nomenological adjustment of the model parameters [p[. 
While fitting the parameters one has to decide for one of 
the current forms of the center-of-mass correction. The 
first rule to be obeyed is then that all later applications 
should employ_precisely that recipe which had been used 
during the fit Q , but the influence of the actual recipe goes 
further than that. The various approximations in them- 
selves do have slightly different trends with mass number. 
These differences can be counterweighted to a certain ex- 
tend by slight readjustments within the model parameters 
relevant for those trends, but the flexibility of the models 
is limited and systematic differences appear for larger ex- 
trapolations, in case of the center-of-mass correction the 
computation of nuclear-matter properties and finite nu- 
clei at large deformation, e.g. in superdeformed states or 
fission. It is the aim of this paper to present a thorough 
investigation of these subtle side effects from the differ- 
ent recipes for the correction for spurious center-of-mass 
motion which will be simply denoted as cm. correction 
in the following. There are, however, several other correc- 
tions for spurious motion and broken symmetries to be 
made, but the cm. correction is among the most impor- 
tant ones since it is present in all nuclei, while the influence 
of, e.g. the rotational and vibrational corrections can be 
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suppressed by choosing spherical nuclei with stiff potential 
energy surfaces. In those, the rotational correction van- 
ishes and the admixture of vibrational excitations to the 
independent-quasiparticle ground state can be assumed to 
be negligible. 

The paper is outlined as follows: In Section we ex- 
plain briefly the underlying mean-field models in this study. 
In Section we summarise the currently used approxima- 
tions for the center-of-mass corrections and in Section ]^ 
we present and discuss typical observables we find to be 
sensitive to the treatment of the correction for center-of- 
mass motion: binding energy systematics and deforma- 
tion energies. An Appendix presents the formulae needed 
to calculate the cm. correction in relativistic and non- 
relativistic models. 



2 Framework 

We investigate the cm. correction in the frameworks of 
the self-consistent Skyrme-Hartree-Fock (SHF) approach 
1^ and the relativistic mean-field model (RMF) 0,^. In 
both, SHF and RMF, models the corrections for spurious 
motion can be treated non-relativistically. 

The numerical procedure represents the coupled SHF 
and RMF equations on a grid in coordinate space using a 
Fourier definition of the derivatives and solves them with 
the damped gradient iteration method |^. We consider 
both spherical and axially symmetric deformed configura- 
tions. 

Pairing correlations are treated in the BCS scheme us- 
ing a delta pairing force |10 l^air — VqS{ri — r2). The 



pairing strengths Vp for protons and Vn for neutrons de- 
pend on the actual mean-field parametrisation. They are 
optimised by fitting for each parametrisation separately 
the pairing gaps from a fourth-order finite-difference for- 
mula of binding energies in isotopic and isotonic chains 
of semi-magic nuclei throughout the chart of nuclei. The 
pairing-active space is chosen to include the number of 
one additional shell of oscillator states above the Fermi 
energy with a smooth Fermi cutoff weight, for details see 

rtnl. 



3 The center-of-mass correction 

The cm. correction - the change in binding energy from 
projection-after-variation in first-order approximation - 
is given by 



^c.m. 



1 



2mA 



iH 



(1) 



We will denote this as microscopic cm. correction through- 
out this paper, referring to the fact that it is calculated 
from the actual many-body state. The explicit depen- 
dence on mass number of p) might cause some trouble 
when calculating mass differences, see iQ. In the frame- 
work of energy density functionals, the factor A should 
be replaced by the integral of the local isoscalar density; 



since we will use (ffl) for correction-after-variation only 
throughout this paper, however, this would make no dif- 
ference. 

Pern = J2k Pfe is the total momentum operator in the 
center-of-mass frame, which is given by the sum of the 
single-particle momentum operators. Although the BCS 
state has vanishing total momentum (Pcm) — 0, it is not 
an eigenstate of Pcnn but has a non- vanishing expectation 
value of its square 



(P?.m.) = II ^'a pL - I] yWp Pa/3 ' P 



•Pa/3 



a,/3 



-I- ^ VaUaVpUp Pa/3 • Pg^g . (2) 

a,/3 

The a and /? denote single-particle states. The p\^ are 
single-particle expectation values of the square of the single- 
particle momentum operator. They appear only in the 
direct term of the correction. The 'Pap are off-diagonal 
single-particle matrix elements of the momentum oper- 
ator. Their squares result from the exchange terms in 
(P^fjj). The further evaluation of (pi) is outlined in the 
Appendix. 

Although S™^ as given by (|l|) is already the first- 
order approximation for the momentum-projected bind- 
ing energy there exist numerous further approximations 
for -E™! in the literature: 

(A) The full correction (g) is considered in the variational 
equation and the calculation of the binding energy _Etot 



5{&n 



^c.m.. 
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^c.m. 



(3) 



This was employed so far only in the fit of the Skyrme 



1J,|1J . For HF states with- 
ditional term to the equa- 



interactions SLy6 and SLy7 
out pairing Ec.-m. gives an ad( 
tions-of-motion but the HF equations can be solved 
as usual. For HFB states the mean field becomes state 
dependent, which requires an additional constraint on 
orthonormal single-particle wave functions in the vari- 
ational equation as described in |lq| . The numerical 
solution of the resulting equations of motion is very 
costly, especially in deformed calculations. Therefore 
this scheme was employed so far only for the descrip- 
tion of doubly-magic nuclei where an HF state can be 
used. 
(B) The cm. correction is omitted in the variational equa- 
tions, but the microscopic correction (||) is considered 
when calculating the total binding energy 
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This a posteriori correction scheme is used, e.g., for 
the Skyrme forces SkIl-SkI5 [^ and the RMF forces 
NL-Z 0, and PL-40 |l9|. 
(C) The cm. correction is approximated by its diagonal 
(direct) terms 
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Table 1. Compilation of nuclear matter properties for a num- 
ber of typical parameter sets. SkM*-SLy6 are Skyrme forces, 
Dl and DIS are Gogny interactions, and NLl, NL-Z, NL3, and 
NL-SH RMF forces, 'scheme' is the scheme for cm. correction 
employed in the fit of the particular interaction, Ovoi denotes 



the volume coefficient or energy per nucleon, 



the (volume) 



symmetry coefficient, and a^urf the surface coefficient. Empir- 
ical values for the volume coefficients derived from the liquid- 



drop model are Ovoi = —16.0 ± 0.2 and Oa; 



32.5 ± 0.5 41 



Force 



Ref. scheme 



Ctvol 

[MeV] 



[MeV] 



f^surf 

[MeV] 



SkM* 


^1 c 


-15.75 


30.0 


17.6 


SkP [ 


221 C 


-15.92 


30.0 


18.0 


SkT6 [ 


|1 C 


-15.94 


30.0 


18.1 


E, 


5 C 


-16.00 


26.5 


18.2 


z. 


51 B 


-15.85 


26.7 


16.9 


Skll [ 


111 B 


-15.93 


37.5 


17.3 


Skis [ 


1(1 B 


-15.96 


34.8 


17.5 


Skl4 [ 


1(1 B 


-15.92 


29.5 


17.3 


SLy4 [ 


yi c 


-15.97 


32.0 


18.2 


SLy6 


141 A 


-15.92 


32.0 


17.4 



Dl 
DIS 



C 

c 



-16.02 
-16.02 



20.3 

18.2 



NLl 


3C 


^ 


-16.42 


43.5 


18.7 


NL-Z 


17 


1 B 


-16.19 


41.7 


17.7 


NL3 


2£ 


1 D 


-16.24 


37.4 


18.5 


NL-SH 


27 


D 


-16.35 


36.1 


19.1 



- where t is the local kinetic density 
before variation 



but employed 
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The contribution from the cm. correction to the (non- 
relativistic) equations of motion has the same structure 
as the kinetic term and leads to a renormahsation of 
the mass of the nucleons 



m TO V A 



(7) 



This scheme is used for most Skyrme interactions like 
SIII 111, SkM* 111, SkP H, the Skyrme forces of 
Tondeur Q, SLy4 and SLy5 ||l|l, the Skyrme inter- 
actions used in semi-microscopic ETFSI calculations 
[|^, and the Gogny forces [|5). 
(D) The microscopic cm. correction (|^) can be evaluated 
analytically for harmonic oscillator states. Using the 
usual parameterisation of the oscillator constant from 
the Nilsson model one obtains an estimate for (nl) as 



770SC 



= -^ Al A-^/^ MeY 



(8) 



This is used for the RMF forces NL-SH [g^, TMl 
[ p8[ , and NL3 |^ . Because (||) does not depend on the 
many-body wave function this gives no contribution to 
the variational equations. 



(E) The cm. correction is approximated by 



E, 



'fit 



-17.2 A' 



MeV 



(9) 



which is a fit to values of the full a posteriori correction 
for the fit nuclei calculated with the Skyrme interac- 
tion Z„ H. This is used for the RMF forces NLl and 



NL2 130 



Some of these approximate schemes for cm. correction are 
discussed and compared in ]31| , where also an additional 
one is proposed which is not used in the models discussed 
here and a short historical overview is given. However, 
the authors of pl| restrict themselves to the quality of 
the various schemes to approximate (nl). 

As can be seen from the harmonic oscillator value (||) 
- which will turn out to give at least the right order of 
magnitude for E™^ in heavy nuclei - the cm. correction 
decreases with increasing mass number A and vanishes for 
infinite homogeneous nuclear matter (which is translation- 
ally invariant). Since the total binding energy increases 
with A, the relative contribution of the cm. correction to 
the total binding energy is largest for very small nuclei. 

In the correction-before-variation schemes (A) and (C) 
the cm. correction directly affects other observables than 
the total binding energy as well, for example single-particle 
energies. Like in case of the total binding energy the effect 
of the correction decreases with increasing A. All other 
schemes deliver an approximate correction of the binding 
energy only, but analogous a posteriori correction of order 
(Pern.) can be performed for the local density p(r) and 
with that all observables derived from it [nj32]. 

We mention in passing that no cm. correction of the 
binding energies is performed in macroscopic-microscopic 
models, but the smooth part of Ec.m. is implicitly con- 
tained in the macroscopic part of these models. 



4 Results and Discussion 

4.1 Nuclear Matter Properties 

From the model systems of infinite nuclear matter (INM) 
and semi-infinite nuclear matter (SINM) at saturation den- 
sity one obtains the leading terms in the nuclear liquid- 
drop mass formula 



E ^ Ovoi A 



J^A + flsurf ^2/3 



(10) 



where A is the mass number and / the relative neutron ex- 
cess given by / = {N — Z)/A. The relation between self- 
consistent mean-field models and the liquid-drop model 
and its refinements like the droplet model is discussed in 
[ P9[ . While the volume coefficient Ovoi and volume symme- 
try coefficient Osym calculated in INM are directly compa- 
rable with the liquid-drop values, the surface coefficient 
flsurf has to be handled more carefully. Usually Osurf is 
extracted from calculations of SINM in semi-classical ap- 
proximation which leads to a value of the surface coeffi- 
cient that is not directly comparable to the liquid-drop 
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Fig. 1. Surface coefficient asm! for a variety of effective inter- 
actions. Lines connect interactions which are similarly fitted. 
Open markers denote forces with a simple scheme for cm. cor- 
rection, while filled markers denote interactions where Bc.m. is 
calculated microscopically within scheme (A) or (B). The inter- 
actions with microscopically calculated -Ec.m. have significantly 
smaller surface coefficients. 



model value [ p9p^ , p^ . Instead the particular Skyrme in- 
teraction SkM* - which is fitted with emphasis on its sur- 
face properties - is often used as a reference when dis- 
cussing the surface coefficient of effective interactions. 

The INM properties of effective interactions have been 
extensively discussed elsewhere, see e.g. [|7|, p^ , p9p0|p[ | and 



references therein. The properties of SINIV! 




have gained a 



42 
50 



44,45 



52 



53 



lot of attention as well, see e.g. 
for non-relativistic models and 
for relativistic interactions. We just want to repeat the 
findings relevant for our study. Table H summarises the 
coefficients of (nol) as calculated from symmetric INM and 
SINM for a number of effective interactions. Most non- 
relativistic (relativistic) interactions agree among each other 
in the values for Ovoi and asym although there are some mi- 
nor differences and a few exceptions, but relativistic and 
non-relativistic interactions differ significantly, the RMF 
forces give Ovoi and Osym larger than the empirical values, 
while the Skyrme forces work much better in that respect. 
It is important to note, however, that not all the values for 
Skyrme forces are predictions. For some parametrisations, 
i.e. SLy4, SLy6, SkP, SkT6 and SkM*, INM properties 
were used as input data during their fit. 

The values for the surface coefficient Osurf differ sig- 
nificantly among all effective interactions, see also Fig- 
ure 1^. We present here values calculated with the extended 
Thomas-Fermi method taken from |53|] or calculated with 
the code of M. Brack |5^. The interactions using i?™m 
have systematically smaller surface coefficients Ogurf than 
most of the other interactions. This hints at a correlation 
between the scheme for cm. correction and the surface 
tension of an effective interaction, which becomes obvious 
when looking at pairs of effective forces which are fitted 
in exactly the same way but use different schemes for the 
cm. correction: the non-relativistic Eg. and Z„ , SLy4 and 
SLy6 and the relativistic NLl and NL-Z. In all cases the 
interaction employing E™^ has a value of agurf smaller by 



roughly 1.0 MeV compared to the interaction employing 
a simpler scheme for the cm. correction. 

Although the Skyrme force SkM* employs the simple 
scheme (C) for the cm. correction, it gives also a rather 
low value for Ogurf • As already mentioned, SkM* is an ex- 
ception among the interactions discussed here because it 
was aimed at giving better flgurf via fitting the fission bar- 
rier of ^^°Pu after it was found that all earlier Skyrme 
interactions give wrong values for the surface coefficient 
[ pip9| |. SkM* still serves as a reference for the surface 
properties of effective interactions for self-consistent mod- 
els. Similar emphasis on proper INM and SINM properties 
was taken for SkP [g^, again resulting in a low value for 
flsurf • AH Skyrme interactions employing E'^^ for the cm. 
correction, however, give a similar (small) value for Osurf 
as SkM* which leads to the conclusion that effective inter- 
actions using the microscopic cm. correction give a rea- 
sonable value for Cgurf without taking explicit information 
on the surface coefficient into account when fitting them 
while interactions with simpler cm. correction need addi- 
tional information on the surface tension to give proper 
surface properties. 

This is confirmed when looking at another pair of sim- 
ilarly fitted interactions: the Gogny forces Dl and DIS. 
Both are fitted with the diagonal approximation i^^ m fo^' 
cm. correction, scheme (C), but DIS is a refit of the orig- 
inal interaction taking again explicit information on the 
deformation energy of heavy nuclei into account pq ] . 

In the remaining part of the paper we want to investi- 
gate the reason for that correlation and give some typical 
examples for the infiuence of the difference in the cm. cor- 
rection and the surface coefficient on observables. To that 
end, we have chosen from all parametrisations discussed 
above a pair of Skyrme interactions, i.e. SLy4 and SLy6, 
and a pair of RMF forces, i.e. NLl and NL-Z, which are 
fitted in the same way but use different levels of approx- 
imation for the cm. correction. The first pair compares 
scheme (C) with scheme (B) while the second pair com- 
pares schemes (E) and (B). As already mentioned, the 
correction before variation used for SLy6 becomes very 
complicated in connection with BCS states. Therefore we 
follow the suggestion of Mm and use the a posteriori mi- 
croscopic correction scheme (B) instead which was already 
done in MM- 



4.2 Binding energy of spherical nuclei 

The fit of effective interactions for self-consistent nuclear 
structure models is usually performed with spherical doubly- 
magic or semi-magic nuclei. The origin of the correlation 
between surface tension and the scheme for cm. correction 
found above can be understood looking at the systemat- 
ics of the cm. correction throughout the chart of nuclei, 
see Figure Q, which shows Ec.m. drawn versus the mass 
number A calculated in spherical symmetry for even-even 
nuclei between the (calculated) two-nucleon drip-lines. 

Looking at the results for NL-Z, one sees that the cm. 
correction shows pronounced shell effects and has maxi- 
mum values at shell closures. E™}!: is a measure for the 
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Fig. 2. Amplitude of the correction for center-of-mass motion 
Ec.in. in spherical calculations with the RMF force NL-Z (up- 
per panel), and the Skyrme interactions SLy6 (middle panel) 
and SLy4 (lower panel). The results are scaled with A~^'^ 
to remove the smooth trend of the microscopic cm. correction 
with A. Solid lines connect nuclei in isotopic chains, the dashed 
line is the fit (|9|) used for NLl, the dotted line the harmonic 
oscillator estimate (M). The diamonds denote the nuclei used 
in the fit of the particular interaction. 



localisation of the many-body wave function and in magic 
nuclei the density distribution is indeed somewhat more 
localised: rms-radii and surface thicknesses are smaller 
than the value obtained from interpolating between ad- 
jacent nuclei, which in turn leads to a larger value for the 
variance of the total momentum (P^ ^^ ) . It is to be noted 
that the results from calculations in spherical symmetry 
plotted in Fig may not be taken too seriously when one 
is far away from proton and neutron shell closures: there 
the nuclei can be expected to be deformed with a more 
localised density distribution and somewhat larger Ec.m.- 
However, this does not affect the conclusions drawn here. 



Table 2. Coefficients fivoi and fisurf of the liquid-drop ex- 
pansion (y_l|) of the difference between -E™^. and S^.m. cal- 
culated for the fit nuclei of some Skyrme forces which use 
the diagonal approximation scheme (C). For those parame- 
terisations where an equally fitted interaction with full micro- 
scopic cm. correction exists, Aa^oi (^"surf) gives the differ- 
ence between the volume (surface) coefficients of those interac- 
tions, e.g. ziaaurf(SLy4) = asurf(SLy4) — asurf(SLy6). Most of 
the Skyrme interactions using scheme (C) which are listed in 
Table m are mainly fitted to nuclear matter properties (and to 
a few finite nuclei only) and therefore omitted. 



Force 


Ovol 

[MeV] 


fisurf 

[MeV] 




AUvol 

[MeV] 


ZiOsurf 

[MeV] 


SkM* 

E, 

SLy4 


-0.14 
-0.13 
-0.15 


1.2 
1.2 
1.2 


Eo-- Zc 
SLy4-SLy6 


-0.15 
-0.05 


1.3 

0.8 



™^ oscillates around E°^^ , 



For heavy nuclei A > 100 i?, 
while E^\^ gives a better approximation to E™^ for small 
nuclei up to ^°Zr, but overestimates i?™^ of heavy nuclei. 
It is noteworthy that none of the algebraic schemes for 
the cm. correction used so far is able to approximate the 
smooth trend of E^^ reasonably for all nuclei throughout 
the chart of nuclei. Improved fits are of course possible, 
but the shell effects visible in Fig S can be described only 
when E^^ is employed for the cm. correction. 

At a first glance the results for SLy6 look very similar 
to those obtained with NL-Z, but there are some signif- 
icant differences: The shell oscillations of E"^^ for small 
nuclei are much smaller in the SHF than in the RMF and 
i?™^ from SLy6 decreases more slowly with increasing 
mass number than E^^ from NL-Z. This explains why 
E^^^^ used for the RMF' force NLl is so far off E^^^l for 
NL-Z: The fit (O) was performed using values of E^^ 
calculated with the Skyrme interaction Z^ which gives re- 
sults very similar to those for SLy6 but which are again 
too large compared to values calculated from RMF model 
wave functions. 

The results for E'^.m. ^^^ed for SLy4 are qualitatively 
and quantitatively very different compared to E^^ used 
with NL-Z and SLy6: the shell fluctuations have nearly 
vanished and - even more important - the diagonal part 
E^^^ does not scale with A~-^/^. This demonstrates that 
the exchange contributions to E™^ are of the same order 
as the direct terms - but with opposite sign - and cancel 
them nearly which gives the A dependence of -E™^ found 
for NL-Z and SLy6. 

The wrong trend of E^ m i^^ compared to -E™^ ) with 
A has to be compensated by the effective interaction fjnt 
in order to obtain the proper binding energy of the (fit) 
nuclei. This is the key to understand the large differences 
in the surface tension found in Sect. 4.1: The error in the 



A dependence of simple schemes for cm. correction puts 
a wrong A dependence into the effective interaction which 
reveals itself in the nuclear matter properties. 



It is instructive to fit the difference AEc.m. = E. 



dir 
cm. 



E'^IS between the two schemes for cm. correction for the 
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Table 3. Coefficients of the liquid-drop expansion (|f3) of the 
difference between i?™^. and the approximation to it used for 
some relativistic forces. To obtain a fit of the same quality as 
in case of the Skyrme forces discussed above a curvature term 
has to be included. 



Force 



flvol 

[MeV] 



Q-surf 

[MeV] 



ttcurv 

[MeV] 



NLl 
NL3 

NL-SH 



-0.12 
-0.11 
-0.10 



1.2 
1.1 
1.0 



-3.2 
-2.7 
-2.4 



nuclei used in the fit of the interactions SkM*, E^, and 
SLy4, with the simple mass formula 



^Er,r,_ — 



Ovol 



A 



rf -^ 



(11) 



The resulting values for the parameters Ovoi and Csurf 
are given in Table g together with the actual differences 
Z\avoi and Aasurf between the nuclear matter properties 
of pairs the similar fitted interactions Eg- and Z^, respec- 
tively SLy4 and SLy6. First of all it is remarkable that 
the coefficients Ovoi and Cgurf are very similar for all in- 
teractions, which means that the amplitude of the cm. 
correction is very similar for all non~relativistic interac- 
tions. In the ideal case the nuclear matter properties of an 
effective interaction fitted with a simple cm. correction 
would be larger by the values of the coefficients of ( pT| ) 
than the nuclear matter properties of a force fitted in the 
same way but using E^^ . This is indeed the case for E^ 
and Zcr where Ovoi ~ ^Ovoi and Csurf ~ ^Osurf- However, 
for other interactions things are complicated by the fact 
that several nuclear matter properties were constrained 
during the fit, hiding (but not curing) the problems caused 
by the simple schemes for cm. correction. The constraint 
on Cvoi for SLy4 and SLy6 diminishes the difference Auvoi 
and probably with that also the difference Z\asurf but one 
can expect that there are instead larger differences than 
necessary in the higher-order terms of the liquid-drop ex- 
pansion since AEc.m. itself is not affected by constraints 
on nuclear matter properties. Including symmetry terms 
in ( |l l| ) neither improves the quality of the fit nor signifi- 
cantly changes the values of avoi and asurf ■ This simple pic- 
ture works well to explain how the use of the approximate 
cm. correction scheme (C) affects nuclear matter prop- 
erties of non-relativistic interactions. The fits of AEc.m. 
with (y_l|) for these interactions give an average deviation 
oUAEl^ « 0.09 Me V^ 

For the RMF, however, the usual approximations are 
schemes (D) and (E). Here, the mass formula ( [ll| ) gives 
very bad fits of AE^.m. with SAEl^ « 0.6 MeV^for all 
interactions. To obtain a fit of the same quality as in the 
case of non-relativistic interactions a curvature term has 
to be included 



AEc.m. — 



flvol 



A- 



^surf 



^2/3 



.^1/3 



(12) 



The volume and surface coefficients are of similar size as 
in case of the SHF model. In the RMF we have a pair 
of similarly fitted forces, namely NLl with the recipe (E) 
and NL-Z with the microscopic cm. correction (B). And 
again, the differences of nuclear matter properties between 
these forces are quite close to the coefficients hi for NLl: 
flvoi ~ ^flvoi = —0.23 MeV and agurf ~ ^Osurf = 1-0 MeV. 
Altogether we find that the RMF behaves much similar as 
the SHF concerning the impact of the cm. correction on 
extrapolation to nuclear matter properties. 

One may have the impression that the correlation be- 
tween the recipe for the cm. correction and the emerg- 
ing surface tension is a merely technical problem, but one 
has to be aware that the choice of the recipe has remark- 
able consequences for extrapolations. Besides the different 
extrapolation of total binding energies visible in Fig. m, 
we discuss in what follows the dramatic consequences for 
strongly deformed systems. 



4.3 Deformation Energy of Heavy Nuclei 

We have already seen that a different treatment of the 
cm. correction in the fit of an effective interaction leads 
to very different surface properties even when the inter- 
actions are otherwise fitted in exactly the same way. To 
get an impression of the amplitude of this effect, we now 
look at superdeformed states and fission barriers of heavy 
nuclei, where experimental information about the defor- 
mation energy is available even for very large deforma- 
tions. As an example we have chosen ^^"Pu which pro- 
vides the standard testing ground for the capability of 
self-consistent nuclear mean-field models to describe fis- 
sion barriers, see e.g. |l|,|6|J5|,|5^,||,|9) . 

In the following we will present deformation energy 
curves calculated in axial symmetry, allowing for reflection 
asymmetric shapes with a damped quadrupole constraint; 
for numerical details see [^QJ. The deformation energy 
curves are shown versus the dimensionless quadrupole mo- 
ment of the mass density which is defined as 



/32 = 



47r 
SArg 



''Y20) with ra = 1.2 A^/^im. 



which yields again an acceptable error of SAE^ 



2 

cm 



and which has to be distinguished from the generating de- 
formation parameter which is used in macroscopic-micro- 
scopic models [ p3[ . 

The upper panel of Fig. shows the deformation en- 
ergy of ^''"Pu calculated in the SHF model with the inter- 
actions SLy4 and SLy6. The first barrier is slightly lowered 
for triaxial shapes which are not considered here, but this 
has no influence on the conclusions drawn in this paper. 
The authors of ^^, however, report a reduction of the 
first barrier of 2.1 MeV for SLy4 when allowing for tri- 
axial shapes. For deformations smaller than the value at 
the isomeric state at (32 — 0.85, the minimum configura- 
tions turn out to be reflection symmetric, while for larger 
deformations the flssion path prefers asymmetric shapes. 



0.09 MeV . Both Skyrme interactions give very similar deforma- 



The actual values of the coefficients di are listed in Table g|. tion energy at small deformations inside the flrst barrier. 
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Fig. 3. Deformation energy (normalised to the prolate ground 
state) of ■^'*°Pu calculated with the Skyrme interactions SLy4 
and SLy6. The lower panel shows the change of the cm. correc- 
tion compared to the ground-state value in the same scale as 
the deformation energy in the upper panel. The various shapes 
along the fission path are indicated by the contours of the den- 
sity at po = 0.07 fm"^. 



but beyond the first barrier a significant difference be- 
tween the deformation energy curves of SLy4 and SLy6 
becomes visible which increases steadily with deforma- 
tion. At the second, superdeformed (SD) minimum the 
difference is already 1.9 MeV and increases to 6.3 MeV 
at /32 = 2.2. SLy4 gives a much broader and higher fission 
barrier than SLy6 which will make a huge difference when 
calculating fission half-lives. The difference between the 
interactions is not caused by the variation of the cm. cor- 
rection with deformation as can be seen in the lower panel 
of Fig. m it is caused by the difference in Cgurf between the 
two interactions. This was already noticed by the authors 
of Ref. |l|. 

The excitation energy of the SD minimum is deter- 
mined mainly by the interplay two very different effects: 
The variation of surface (and surface symmetry) energy 
with deformation contained in the bulk properties of the 
effective interaction (the Coulomb energy which varies 
with deformation as well is nearly the same for all interac- 
tions) and the variation of shell effects with deformation 
(which is fixed by the shell structure at spherical shape). 
SLy4 and SLy6 give nearly the same single-particle ener- 
gies for spherical nuclei, the fission paths are identical, the 
difference visible in Fig. |^ is caused only by the difference 
in surface tension between the two interactions. 

This gives a possible explanation for the finding of the 
authors of |Q that the Skyrme interaction SLy4 overes- 
timates the excitation energy of the SD minima in nuclei 
around A Ki 190 and A « 240 while it reproduces nicely 
the separation energies within the first and SD minima. 
In each region either all ground-states and all SD minima 
have roughly the same deformation, therefore the sepa- 
ration energy calculated between states in the same well 




Fig. 4. Fission barrier of Pu calculated with the RMF 
forces NLl and NL-Z drawn in the same manner and using 
the same scales as in Fig. H. 



is not affected by the surface tension. This might explain 
also to some extent the preference for SkM*, SkI3 and 
SkP for the description of superdeformed states in ^^^Hg 
found in p%. All three forces have rather small surface 
coefficients asmf (see Table |^ and Fig. |l|), SkI3 because it 
employs the microscopic cm. correction, SkM* and SkP 
because the value of Osurf was constrained during the fit 
of these particular interactions. However, the excitation 
energy of superdeformed states is determined by the in- 
terplay of surface properties and shell structure, both have 
to be properly described to reproduce experimental data 
throughout the whole chart of nuclei. More investigation 
in that direction is needed. 

Figure ^ shows the same as Fig. ^, but for the RMF 
forces NLl and NL-Z. The qualitative features of the fis- 
sion barrier are the same as those found for the Skyrme 
forces, but there are significant quantitative differences. 
Comparing both figures it can be clearly seen that all 
barriers are much smaller for the RMF forces compared to 
the Skyrme interactions. Again, the RMF force which em- 
ploys the "cheaper" cm. correction gives the larger defor- 
mation energy compared to the ground-state value. There 
are now some differences between NLl and NL-Z at small 
deformations as well, but for large prolate deformations 
we regain our finding for the Skyrme forces: The differ- 
ence in the deformation energy increases with deforma- 
tion, it is already 1.9 MeV at the SD minimum and 4.7 
MeV around (32 = 2.2. The RMF forces shift the SD min- 
imum to a slightly larger deformation f32 ~ 1.0 than the 
SLycc forces. For NL-Z the excitation energy of the SD 
minimum of —0.8 MeV is definitively too small compared 
to the experimental value of 2.4 MeV, while the value of 
2.8 MeV predicted by NLl is quite close. We expect, how- 
ever, that the better description of the SD state obtained 
with NLl is accidental because the surface coefficient of 
this interaction is definitively too large. 
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Again the difference of the deformation energy curves 
is not caused by the actual contribution of the cm. cor- 
rection to the binding energy, see the lower panel of Fig. ^. 
E^^ varies only by a few 100 keV along the fission path 
while the value for E^\^ is of course constant. 

We have to revise our conclusions from p3| where the 
fission barriers of ^"^^Pu calculated with PL-40 |I9| and 
NLl were compared. Both interactions are fitted in the 
same way but differ in two details: the treatment of the 
cm. correction - where PL~40 employs the same scheme 
as NL-Z - and the density dependence, where NLl and 
NL-Z employ the standard non-linear self-interaction of 
the scalar field while PL-40 uses a stabilised self-interaction, 
see [0,@ for details. We have concluded in Ref. |Q that 
the difference between the fission barriers is caused by the 
difference in the density dependence of the two forces PL- 
40 and NLl, but as can be seen comparing Fig ^ with the 
results of [q9| NL-Z and PL-40 give very similar fission 
barriers, therefore the difference between NLl and PL-40 
is caused by the different treatment of the cm. correc- 
tion during the fit of these interactions. This means in 
turn that using the stabilised non-linearity instead of the 
standard non-linearity in the RMF has only a very small 
influence on the fission barriers of heavy nuclei. 

It is to be noted that the mere cm. correction (^ 
becomes inconsistent in the asymptotics of fission (i.e. be- 
yond the scission point). One still would use one common 
-E-c.m. whereas the two fragments should now acquire each 
their own i?c.m.- The problem can be resolved by includ- 
ing also the correction for spurious quadrupole vibrations 
p0| , |61[ . But this complication can be neglected when look- 
ing at the fission barrier of heavy nuclei where the scission 
point is far beyond the fission point. 



4.4 Deformation Energy of Light Nuclei 

Figure ^ shows the deformation energy for the nucleus ^"'S 
calculated again with SLy4, SLy6, NLl, and NL-Z. The 
potential landscapes for the relativistic and non-relativistic 
interactions look very different: SLy4 and SLy6 predict a 
rather soft potential energy surface with a shallow oblate 
minimum, while NLl and NL-Z predict an energy surface 
with deep prolate and oblate minima of comparable bind- 
ing energy (which might be connected by triaxial configu- 
rations). The main difference between the models appears 
at small deformations around /32 = 0, while all interac- 
tions give similar predictions for the deformation of the 
prolate and oblate minima respectively. This is interest- 
ing for itself as it has dramatic influence on the excitation 
spectrum of this nucleus, but we do not want to discuss 
this further here and refer to pg and references therein. 
In contrast to the heavy system ^"^^Pu there now ap- 
pear visible differences at small deformations as well, which 
are directly related to the different variation of E^-^ and 
-^c"m with deformation, see the lower panel of Fig. pi Com- 
paring SLy4 and SLy6, one sees that the diagonal part of 
the cm. correction oscillates with larger amplitude than 
the full correction including exchange terms. The relativis- 
tic forces NLl and NL-Z show a small difference at small 




Fig. 5. Deformation energy (upper panel) and variation of 
the cm. correction (lower panel) of *'*S calculated with NLl, 
NL-Z, SLy4 and SLy6. The deformation energy E and the cm. 
correction Ec.m. are normalised to their values at the prolate 
minimum. The two panels use again the same energy scale. 



deformations as well, but this difference has another ori- 
gin: -E™^ is slowly varying while the fit value is simply 
constant. At large prolate deformations the difference in 
surface tension between SLy4 and SLy6 gives again a large 
difference between the potential energy surfaces. The dif- 
ferences caused by the cm. correction among these par- 
ticular relativistic and non-relativistic forces are rather 
small compared to the difference between the predictions 
of the non-relativistic and relativistic models and it is to 
be noted that SLy4 and SLy6 are by no means represen- 
tative for the predictions of Skyrme interactions, see [ p8[ . 
However, the differences we see for small nuclei at small 
deformations caused by differences in the variation of the 
cm. correction with deformation are subtle compared to 
the huge effect caused by the difference in surface tension 
seen at large deformations. 



5 Summary 

We have discussed the various approximate schemes for 
performing the cm. correction in connection with self- 
consistent mean-field models. In particular, we have scru- 
tinised the effect of the actual recipe for the cm. cor- 
rection used while fitting a parametrisation. The basic 
point is that the different recipes differ in their trends 
with mass number A and isospin. Mismatches thus inher- 
ent in a recipe are counterweighted to a certain extend 
by (automatic) readjustments of the model parameters. 
This, in turn, builds into the forces different properties 
concerning extrapolations. 

We have considered the cm. correction of order (P^ ^ ) 
throughout. The full evaluation of -E™^ was compared 
with the approximation through the diagonal part of that 
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expression and with a simple estimate oc A''^ . The com- 
parison of the cm. correction as such shows significant dif- 
ferences between these three schemes. The full correction 
shows considerable shell effects which are completely ab- 
sent in all simple estimates. Moreover, the average trends 
of Ec.m. with A and I of these three recipes differ, which 
has consequences for the properties of effective interac- 
tions fitted with simple schemes for cm. correction. We 
see a strong effect on the surface energy where the forces 
with approximate cm. correction (diagonal part or simple 
estimate) give systematically larger values than interac- 
tions employing E^^ . 

This error in surface energy leads to a much different 
evolution of energies with deformation. As a consequence, 
fission barriers in actinide nuclei are about 4 MeV larger 
for those forces which employed the approximate recipes 
during their fit. Furthermore, predictions of the relative 
height of deformation isomers and collective vibrations 
therein will come out significantly different. For light nu- 
clei - where the relative importance of the cm. correction 
is larger - one even has to be aware of differences at small 
deformations. Altogether, we strongly recommend to use 
consistently the full correction (P^ ,-„ )/(2toA) and forces 
derived with that recipe. Any further approximation be- 
yond that (already approximate) stage induces systematic 
errors in extrapolation to heavy nuclei, large deformation, 
and probably also to nuclei far off the stability line. 
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A The Calculation of (P^^ ) 



In this appendix we present the formulae needed to cal- 
culate the expectation value of (Pc.m.) i'^ relativistic and 
non-relativistic models assuming spherical, axial symme- 
try or triaxial symmetry for a BCS many-body state. Ex- 
pressing (||) in terms of matrix elements Aaa of the Lapla- 
cian and Vq.^ of the nabla operator respectively, one ob- 



tains 



\-^ cm./ 



-VaUaVpUpVafS-'^a^)] ■ (13) 

In case of a time-even many-body state this simplifies 
with Aaa = Aaa, Vg^ = V*^ and Va/3 = -'^ap to 

(P?.m.) = -2^1^ I ^ vl Aaa + ^ ^aVf3{VaVi3 + UaU/s) 
Q>0 a,/3>0 

2^ 



x(|V„0p + |V„^PjJ . (14) 

A.l Non-Relativistic Models 

A. 1.1 Spherical Symmetry 

In spherical symmetry non-relativistic spinors are given 

by 

^ajim{r) = {r\aj£m) = tpajeir) (ije,n{^, 4>) (15) 

where a is the principal quantum number, j the total an- 
gular momentum, I the orbital angular momentum, and 
m the projection of the total angular momentum, ij^aji is 
the radial wave function to be treated numerically, while 
the spinor spherical harmonics fijim pq | can be treated 
analytically. When calculating the cm. correction (|lj) 

l2 



\-^ cm./ 



+J 



fc2 Y^ 



y {aj£m\A\ajim) 



m=-] 



" ^_^ '"ajt- '"PJL (Vajl VpjL + Uajg Ujjjl) 
a,f3,j,J,£,L 



xE 



E 



{ajem\W\(3JLAI)\ 



(16) 



the summation over m and M can be performed analyt- 
ically: the sum over m (for given j and £) of the matrix 
elements of the Laplacian gives simply 2j + 1 times the ex- 
pectation value of the Laplacian of ^paji, while the square 
of the matrix elements of the nabla operator reads 



+j +J 

E E 

m——j M—~J 



\{aj£m\y\l3JLM)\ 



= (-1)^+^(2 J +l)(2j + l) 



J ] 1 



X [L 5l.L-l Aap Bfja + i (5l,£-1 Ai3a ^a/s] (17) 



ith 



A 



Q/3 



dr r' 



Baf3 ^ I dr 






(18) 
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The square of the 6j symbol occurring in ( |17| ) is given by with 



^L^j (2j + l)(2J+l)(2£ + l)(2L+l) ^ I ^^" ^'^ 



a. 11 



x((5j/+i Sjx+i + 5j,e-i Sj,L-i) 



(J_^+3)(^_j+3) ^^0- i^.^u,,^^ v-^^ 



A^-p^ jd\Y.a^^i--^^V^^^ 



(2j + l)(2J+l)(2^+l)(2L + l) 



x(^,,+r 5 J,,., + ^,,_r ^.,.H-i) (19) ^"/^ = /d^'^E '/ V^i'^^^ V <") 



a.Tj 



A.1.2 Axial Symmetry ^«^ = /dV^ a ^(— ^') V <") (27) 

CT,7? 



A. 2 Relativistic Models 



Assuming axial symmetry, a spinor with angular monren- 
tum projection m is given by 

\Wa {'TjZ) exp \i\7n+ 2JV\ J Relativistic kinematics plays no role for nuclear ground 

states, the cm. correction can be calculated in non-relativistic 
The V'a can be chosen to be real. The matrix elements approximation in relativistic models. The main difference 
of the nabla operator needed to calculate (|l|) or (|l|) are between relativistic and non-relativistic models is that the 
given by relativistic spinors have four components which changes 

V T7* -A A A the calculation of (P2^.}. 

\7 Y7 — A A A A. 2.1 Spherical Symmetry 

-T- 2\""ic.mf,-i ^ap ^ap ~^ ""^"^^-0+^ ^aP ^apj The rclativistic spinors in spherical symmetry are given 
Vo/5 • V:^ = i <5,„„, 1 <5„^, 1 C„;3 C„^ (21) by 

with / ihith 



^a/3 = y'dV^^i-)9,4'^) (22 



<?_(r) = (r|a.r.) = . fc/;);^-™^;^ \ (28) 

iV'QK'(r)l2_K,„(?9, .^) y 



The quantum number k is related to the total and orbital 

B^p = (d?r Y^ Vi"^ (dr ± ""^ ^ M V^*"^ (23) angular momentum by j = |k| - i and ^ = j + i,^. 

•^ cr=± \ ^ The real functions ip'^'^'^ are the radial part of the upper 

/.^ / jYL„ j^ 2.\ ._. and lower components of the wave functions, the i^^.m 

d '^ ^^ {~'^)V'a I C'r H - I "0/3 (24) are spinor spherical harmonics |Q with angular momenta 

<^=^ corresponding to k. The variance of the total momentum 

1 ,1 , 1 , ■ • ^ ^-i ^11 is now given by 

where the volume element is given by d r = Znarraz. 

+j 

A 1 o /- 4- ■ D 4. 4.- (Pr m ) — ^^^ > v'i^. > (aK?Ti|Z\ laKm) 

A. 1.3 Cartesian Representation ^ ^■™-' /-^ "'*■ Z-^ ^ ' ' ' 

a,K m——j 

In Cartesian representation the non-relativistic spinors _;j2 y-^ ^^^ ^^^ ^^^^ ^^^ _l_ ^^^ ^^^^ 

are given by ^_^ ,^^ 

'Z'.(r)=f4!;^Jj'"J+^^^!:Ij'"j') ' (25) xE i: |(a..m|V|/?XM)|^ (29) 

VV^" (O+lV-a HI-)/ m=-jM=-J 

where the ipa are real functions. The matrix elements with 

of the nabla operator needed for the calculation of ( [13|) 

are given by +■' /■ 

2 2 E (««"i|^l««"i) = (2j + 1) 5] / dr r2 V^(:j) Z\ ^^ 

V./5 • K^ - K-, + K^ (26) (30) 
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and 



+3 

>; 


+ J 

> ; {ajK\y\pKM) ' ( 

M=-J 




X 


L5^M^i^M_i A^l^ Bf^^ + ^^i(„)_f(r,)_i A^^ B^p 



where A and B are given by 






dr ' 



V'q« ( (^r 



L'''' + 1 



^x 






dr r^ ^„, (9, - i^) ^/3^ . (32) 

The 6j symbol appearing in (p^) is again given by (19). 



A. 2. 2 Axial Symmetry 

The relativistic spinors in axial symmetry are given by 



<?a(r) 



( i^i++^(r,z)exp[i(m-i)0]\ 
ip'i^~\r,z) exp [i(m+ \)(f\ 
Wa {r, z) exp [i(m - i)0] 



(33) 



\i-0i ^(r, z) exp [i(m+ 5)0] / 



where (r, z, 0) are cylindrical coordinates and m is the 
projection of the total angular momentum. The i/;''''^' are 
real functions, where 77 denotes upper and lower compo- 
nents and cr/2 the spin projection. The matrix elements 
entering (21) are now given by 



^aoi 






dt 



■dr- 



(m^^,; 



«' 



{ria) 



A 



a/3 



dV^V^^)9.^('"^) 



7^,(T 



^a/3= /dV^V^i"-) 9.± 



m/3 



Co/; 



77 cr 



^r^ 



mp 



V 



(';,■ 



(34) 



A. 2. 3 Cartesian Representation 

The relativistic spinors in Cartesian representation are 
given by 



/Vj 



(+++) 



<?o(r) 



(r) +1^;, 



(++-) 



(r)\ 



^i+-+^(r)+i^i+— )(r) 
V^++^(r)+i^i-+- 
Vvi""+^(r)+i^i ^(r)/ 



(35) 



where the ipa^ are real functions. 77 and a have the same 
meaning as in the axial case, g denotes real and imaginary 
part of a spinor component. The matrix elements of the 
Laplacian are given by 



A^^ = Jd\Y.iji^^s)^^in^e) 



(36) 



while the integrals needed to calculate the matrix elements 
of the nabla (pq) are given by 



Aq/3 — 



dV^ Vi^'^^V^""'^) 



ry.V.Q 



dWe' 



{v,-o-,s) ■ 



<y,ri,Q 



0/3- /'dV^V^'^'-'^^v 



(■nag) 



B 



4"-^) 



■qag) 



(37) 



iy,n,e 
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